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Abstract 


This  paper  develops  optimal  statistical  estimation  formulation  for  multi  - 
dither  adaptive  optics  control  loops  which  potentially  can  enhance  stable  beam 
control  performance  in  the  presence  of  spurious  signal  like  noise.  A multi- 
dither  autofocus  example  is  chosen  to  compare  estimated  assisted  and  conven- 
tional adaptive  optical  performance  in  the  presence  of  speckle  generated  noise. 


I.  Introduction 


Adaptive  optics  has  truly  become  an  important  subficld.  As  the  ideas  in 
adaptive  optics  cross  interdisciplinary  lines  of  both  the  optics  community  ;uid 
the  control  community  (although  perhaps  not  universally  recognized),  it  is 
imperative  that  new  techniques  from  each  community  evolve  to  enhance  the 
pi  r fomance  of  adaptive  optical  systems.  This  paper  represents  one  of  these 
iuderdisoipl inary  crossings  in  order  to  utilize  experience  and  techniques  from 
the  control  community  to  enhance  the  field  of  adaptive  optics. 

Adaptive  optical  systems  are  inherently  closed  loop  systems  whereby  a 
feedback  signal  is  used  to  detect,  system  error  in  some  manner . This  error 
detection  is  then  used  to  control  an  adaptive  optic  in  order  to  enhance  overall 
optical  performance.  The  inherent  performance  of  closed  loop  control  lies  in  its 


ability  to  correct  for  a process  that  may  not  be  as  accurately  determined  as 
one  would  like.  However,  the  limitation  of  feedback  control  is  that  the  feed- 
back signals  must  be  as  precisely  determined  as  the  accuracy  level  required  of 
the  controller  [1|.  Since  in  adaptive  optics  the  signal  may  he  corrupted  by 
several  noise  effects  such  as  speckle,  photon  noise,  and  detector  noise,  it  is 
imperative  to  obtain  a filtering  of  this  noise  in  order  to  properly  achieve  the 
potential  system  accuracy.  Without  fully  realizing  the  use  of  optimal  tech- 
niques, the  optics  community  has  actually  benefited  from  a form  of  filtering 
that  is  inherently  present  in  existing  adaptive  optics  servo  hardware.  The 
inherent  rum  log  response  of  these  servos  to  feedback  signals  yield  a filtering 


of  the  signal.  However,  the  filtering  is  at  best  subopt imal  as  it  does  not 


include  the  prior  statistical  knowledge  of  the  phenomenon  perturbing  the  feed 
back  signal.  Nor  does  it  use  any  of  the  available  information  in  an  optimal 
sense.  The  filtering  is  strictly  a bandpass  operation.!  ' y 


l 


'Ihis  paper  develops  the  use  of  estimation  techniques  as  well  known  in  the 
control  community  [2,3,4]  in  order  to  obtain  an  optimum  estimate  of  the  feed- 
back signal  for  multidithcr  adaptive  optics  (see  [5,6]  Coherent  Optical  Adap- 
tive Techniques  and  Image  Comnensation) . The  feedback  signal  is  then  used  in 
a stochastic  approximation  scheme  in  order  to  obtain  a convergent  feedback 
controller  for  the  multidither  adaptive  optics.  The  estimator  will  give  an 
optimal  estimate  of  the  feedback  signal  in  the  presence  of  both  speckle  noise 
and  detector  noise.  The  results  of  the  technioue  are  applied  to  an  auto- 
focusing  scheme. 

The  control  is  based  on  optimizing  the  intensity  as  a function  of  the 
distance  between  the  secondary  and  primary  mirrors  for  a Cassegrain  tele- 
scope. Maximization  logic  is  based  on  the  parameter  gradients  of  the  inten- 
sity. Second  order  gradients  are  not  used  in  this  simple  maximization  logic. 
However,  they  can  easily  be  incorporated  in  multi-element  adaptive  optic 
systems . 

Tire  worse  case  speckle  problem  lias  been  addressed.  This  is  the  additive 
speckle  problem.  In  the  example,  a simple  phase  lock  loop  was  modeled  for 
comparison  purposes . An  ACC  compensation  was  not  simulated  as  it  does  not 
add  any  utility  to  the  models.  It  is  shown  that  the  technique  in  this  paper 
yields  significant  performance  improvement  even  with  ideal  modeling  of  the 
comparison  loop. 

Section  I gives  the  problem  statement  including  a discussion  of  the 
sinusoidal  perturbation  adaptive  controller  and  the  basic  estimation  equations. 
Section  II  considers  the  problem  of  speckle  interaction  on  adaptive  optics  from 
the  standpoint  of  induced  signal  spectrum  components,  and  the  modeling  of  an 
estimator  addressing  this  problem.  Section  111  gives  the  estimator  equations 
in  the  continuous  form  for  possible  analog  implementation.  Section  IV  gives 
the  control  philosophy  in  terms  of  stochastic  approximation.  Section  V demon- 
strates the  results  of  applying  estimation  theory  to  the  autofocusing  problem 
and  Section  VI  contains  the  conclusion. 


II.  Problem  Statement 


The  philosophy  of  multidither  adaptive  optics  is  to  obtain  a maximum  of 
intensity  by  using  a feedback  signal  derived  from  multidithering.  The  intensity 
is  a function  of  adjustable  parameters,  P,  where  P is  a q-vector,  i.e., 

I = UP).  (1) 

The  parameters  are  perturbed  from  a reference  value,  P,  by  the  addition  of  a 
sinusoidal  perturbation  signal,  i.e., 


P = P + 
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The  actual  intensity  at  the  receiver  will  then  be  a function  of  the  perturbations, 
i.e.  , 

i (P)  = up  + n)  ( .U 
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He i . (.'•>)  may  be  expanded  in  a Taylor  series  about  small  amplitude  perturbations 
l|  |fi|  |«|  |P|  ( where  ||(-)||  denotes  the  Euclidean  norm  of  the  vector  (•)  ). 
Tins  yields 
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UP)  = I(P)  + ft  + \ ftT  |pT  n + H.o.T. 

where  H.O.T.  denotes  the  higher  order  terms,  the  superscript  "T"  denotes  the 
vector  transpose, 
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and  where 


3? 1/3P2 j , 3*I/3P23P,,  ....  32I/3P  3P, 


3?I/3P13Pq,  32I/3P2q 


(7) 


Pj  denotes  the  i-th  element  of  the  parameter  vector.  For  small  perturbations 
(in  terms  of  the  norm  of  ft)  the  higher  order  terms  may  be  dropped.  A more 
refined  approach  will  he  considered  when  the  second  derivative  is  employed  foi 
convergence  behavior  monitoring  in  this  paper.  For  the  purpose  of  this 
development,  however,  analog  filtering  will  be  assumed  to  lx?  available  for 
eliminating  the  double  frequency  terms  as  well  as  higher  order  terms.  There- 
fore, if  the  signal  return  were  perfect,  we  would  have  as  the  signal  after 
bandpass  filtering 


up)  = up)  + jp-  n.  (R) 

The  sign  u .1  magnitude  of  the  gradient  of  the  intensity  with  respect  to  the 
adjustable  parameters  are  contained  explicitly  in  the  term  containing  the 
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basic  dither  frequencies.  This  yields  the  information  necessarv  to  make  a 
correction  to  the  adjustable  parameters  such  that  the  intensity  can  be 
increased.  In  particular,  the  new  settings  for  the  parameters  may  be  found 
by  appropriately  choosing  the  desired  intensity  change  at  the  current  time 
and  solving  for  the  required  change  in  the  parameters  by  a steepest  ascent 
algorithm,  i.e., 


Poldk  + f'k  9PT 


wherc  >k  > 0 at  each  time  instant,  k,  and  serves  as  an  effective  gain  parameter 
that  is  controlled  to  assure  reasonable  convergence  towards  the  optimum.  The 
predictor  change  in  I may  be  found  as 


This  is  continued  until 


iM.-  rp  - p ) = c - - — 
DP  1 new  oldj  Lk  3P  3P 
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where  6 is  chosen  for  convergence  control. 

Tlte  above  algorithm  implicitly  used  in  multidither  Adaptive  Optic',  lias 
the  problem  of  the  gradient  so  obtained  is  assumed  to  be  a deterministic 
quantity  and  as  such  is  determined  by  a noiseless  measurement.  Neither  of 
these  cases  are  actually  realized  in  ;tny  application.  It • particular,  in  real 
situations  the  amplitude  of  the  gradient  is  amplitude  modulated  by  speckle  ,:s 
well  as  corrupted  by  signal  like  noise  components,  as  will  lie  showed  in 
Section  II.  Furthermore,  the  measurements  are  not  noiseless,  thus,  the 
returned  signal  from  the  target  may  be  written  as 


v(n  = k(z)U(P)  + 


+ S + S + 


+s  +T+gl  + n 
q * 


where  k(z)  is  the  detector  gain,  z is  the  target  range,  ^1/3P  is  the  gradient 
modulated  at  the  dither  frequency,  s^  is  the  speckle  induced  spectrum  component 
occurring  at  the  i-th  dither  frequency  hand,  T is  the  target  diffusely  distrib- 
uted return,  g is  a random  glint  return  component  and  n is  the  detector  noise. 
The  speckle  induced  spectrum  and  modulation  at  the  dither  frequencies  cause  a 
random  fluctuation  of  the  signal  at  each  dither  band. 

The  second  order  terms  may  be  retained  with  some  further  modeling.  The 
intensity  I ( P)  is  assumed  to  be  representative  of  the  returned  intensity  with  g 
and  T being  zero -mean  processes  for  the  temporal  variation  of  the  glint  and 
diffuse  returns.  In  this  way,  the  temporal  spectral  spreading  of  the  returned 
intensity  may  be  modeled.  Ihat  is,  g and  T are  perturbations  of  the  returned 
intensity  due  to  glint  and  diffuse  changing. 

Thus,  the  measurement  equation  is  as  given  in  Pq . (12).  Since  speckle 
induced  components  which  effect  system  performance  appear  at  the  dither  fre- 
quencies, they  are  indistinguishable  from  the  gradient  vector  components  and 
cannot  be  separated  out  bv  conventional  filtering  techniques.  However,  by  the 
use  of  optimal  estimation  techniques,  it  can  be  discriminated  as  undesirable 
signal  components.  In  fact,  an  estimator  is  structured  in  Section  III  in 
order  to  estimate  the  gradient,  the  speckle  components,  the  glint  component, 
and  the  target  diffuse  component. 

The  structure  of  the  estimator  to  be  used  is  given  without  proof  as  it  is 
well  known  to  the  control  and  estimation  community.  However,  the  structure  is 
given  for  completeness.  The  first  formulation  is  that  of  an  analog  implementa- 
tion of  the  estimator.  The  process  x,  which  is  an  n-vector,  is  assumed  to 
evolve  according  to  t lie  set  of  formal  first  order  vector  differential  cqua- 
t i ons , 

x (t ) = f { x ( t ) } + r,(t)u(t)  (13) 

where  u(t)  is  a zero  mean  white  noise  m-vector  with  covariance 


hlul t)u(i ) ! = Q(t)f (t-i ) 

where  6 ( - 1 ■ a the  Dirac  delta  function.  (The  authors  are  well  aware  of  the  lto 
interpretations  as  well  as  the  general  structure  of  the  result.  However,  it  is 
formalized  for  the  sake  of  the  optics  community.)  lhe  measurement  equation  lor 
the  process  x(t)  is  given  as 


r ■ 1 


I 
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y(t)  = hft)  'x(t)  + v(t) 

where  y may  be  a vector  in  general  but  will  be  assumed  to  be  a scalar  since  the 
measurement  of  the  pi'oeess  of  concern  is  a scalar.  The  vector  v ( t ) is  a zero 
mean  white  noise  process  with  covariance 

b{v(t)v(T)l  = r(t)6(t-T)  . 

It  is  assumed  that  u and  v are  uncorrelated.  The  approximate  conditional  mean 
estimator  for  the  nonlinear  system  is  given  as 


x(t)  = F(t)x  + K(t){y(t)  - hT(t)x(t)) 

A 

where  x is  the  approximate  conditional  mean  estimate, 
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K(tJ 


= P(t)h(t) 
r(t) 
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with 


P(t)  = F(t)P(t)  + Pft)P(t)1  + f>(t)Q(t)C!(t) 1 

P(t)hr(t)h(tJP( t± 
r(t ) 


I 18) 


Die  matrix  P(t)  is  the  approximate  covariance  matrix  for  the  estimation  error, 
x - x,  and  the  last  term  in  l.q,  (18)  represents  the  effect  of  the  measurements 
on  the  differential  equation  for  P(t).  The  structure  is  called  the  extended 
Kalmar  c.  timator  [7,3] . N'ote  that  the  matrices  F,  (1,  Q,  and  r do  not  have  to 
he  time  varying.  It  is  important  to  realize  that  the  estimator  weighs  the 
process  dynamics  and  the  measurements  according  to  their  statistical  knowledge. 
Also,  as  tt,  number  of  measurements  processed  increases,  if  the  estimator  is 
stable,  then  the  knowledge  gained  will  increase.  This  is  significant  since 
some  of  the  signal  processing  algorithms  m use  in  currently  configured  optical 
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systems  lose  information  by  not  taking  this  fact  into  account.  The  discrete 

estimation  algorithm  may  be  easily  established  and  is  not  given. 

The  above  is  an  elementary  discussion  of  estimation  theory.  Tor  more 

details  and  advanced  techniques,  it  is  suggested  that  the  reader  use  the  pre- 

* 

viously  cited  references  for  more  information. 

It  suffices  to  say  that  in  order  to  obtain  an  estimator  that  will 
optimally  estimate  the  quantities  in  Eq.  (12),  it  is  necessary  to  obtain  a 
dynamic  model  structure  as  in  Eqs.  (15-18).  Thus,  dynamic  statistical  equa- 
tions for  each  term  in  Eq.  (12)  must  be  obtained.  This  is  considered  in  the 
next  section. 


III.  Speckle  Effects  Modeling 

/\s  was  indicated  in  the  last  section,  it  is  necessary  to  take  into  account 
all  the  statistical  information  available  about  the  process  to  be  estimated. 

Part  of  the  statistical  informat  ion  is  derived  from  a priori  models  for  the 
time  varying  random  phenomenon.  This  a priori  statistical  information  is 
contained  in  the  models  of  the  form  as  given  in  Eq.  (14).  Thus , in  order  to 
use  the  structure  appearing  in  the  optimal  estimator,  we  must  develop  models 
such  as  these.  This  will  be  accomplished  in  this  section  for  speckle  inter- 
active effects  on  the  adaptive  optics  case  in  question. 

We  are  interested  in  obtaining  a general  expression  for  the  temporal 
power  spectra]  density  (PSD)  for  frequency  components  generated  by  target 
return  dynamic  speckle.  For  the  purposes  of  this  development,  target  back- 
scatter  can  simply  be  characterized  as  having  a specular  and  a diffuse  part. 

Tlte  specular  return  preserves  the  spatial  coherence  of  the  incident  illumina- 
tion. It  cminates  from  target  regions  which  are  normal  to  the  receive!  line- 
of-sight  and  have  a smooth  surface  texture  relative  to  the  dimensions  of  the 
incident  wavelength.  The  diffuse  part  eminate  from  spatially  distributed 
target  regions  within  the  beam  illumination  and  are  characterized  by  the  general 
surface  roughness  properties  and  target  geometry.  Both  Goldfischer  [7j  and 
Crane  [8j  develop  relationships  which  give  the  spatial  PSD,  S(ui,fi)  > in  terms 
of  the  convolution  of  the  object  illumination  function  p(u,v).  We  will  use 
these  relations,  addressing  the  projected  illumination  function  (of  an  other- 
wise three  dimensional  target  shape,  and  resulting  target  illuminated  structure) 
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on 

S (w,fi)  ~ fj’ dudpfu,v)p  (u  - v - ■ ) • (19) 
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in  the  (u,v)  plane,  denoted  p(u,v),  as  the  source  of  both  specular  and  diffuse 
returns.  The  object  shape  dictated  normal  incident  point  is  located  in  Fig.  1 
at  the  (u,v)  plane  origin  for  convenience.  The  remaining  distributed  regions 
in  p(u,v)  account  for  the  object  diffuse  backscatter.  Employing  Goldfischer ' s 
relatively  simple  formulation  for  the  Fraunhofer  case,  we  have 


Or  simply  stated,  the  spatial  PSD  is  the  self  convolution  of  the  target  spatial 
irradiance  (refer  to  |7|  and  |8|  for  proof).  The  spatially  distributed  ampli- 
tude of  p(u,v)  is  directly  relatable  to  the  target  specular  regions  and  the 
overall  backscatter  diffuse  nature  of  the  reflecting  surface.  The  latter  is 
generally  characterized  by  the  fall  off  of  backscattered  intensity  as  a 
function  of  incidence  angle.  The  results  are  a combined  effect  of  the  target 
surface  scattering  properties  and  target  illuminated  shape.  Surface  scattering 
c:u i be  measured  in  the  laboratory  as  Mono- Directional  or  Bi-Directional 
Reflectance  data  (MDR  or  BDR)  for  flatplate  samples  of  the  target  surface 
material,  giving  scattering  behavior  which  is  free  of  shape  effects  (no 
geometrical  effects).  This  point  is  made  to  emphasize  that  the  formulation 
of  p(u,v)  must  contain  the  full  phenomenon  of  surface  and  shape  scattering  for 
any  particular  object  given  a specified  illumination  function  on  that  object. 
Ibis  is  a formidable  task  from  an  analytical  approach  for  complex  target 
structures,  but  reasonably  performed  via  computer  target  shape  decomposition 
techniques  where  laboratory  BDR  data  can  be  used  directly  for  scattering 
calculations.  The  following  discussion  of  shape  responsive  spatial  PSD  uses 
this  object  shape  decomposition  approach.  Ibis  involves  separating  a complex 
structured  target  into  its  fundamental  geometrical  shapes,  identifying  the 
projected  silhouette  of  each  shape  onto  the  (u,v)  plane,  weighting  this 
projection  by  the  beam  illumination  outline,  and  noting  normal  incident 
point  (referenced  as  riie  source  of  target  scattered  glints)  locations  on 
each  beam  weighted  shape  silhouette. 

Ihe  mechanics  of  the  computer  shape  decomposition  approach  for  spatial 
PSD  calculations,  simply  utilize  the  above  stated  relationship  ol  p(u,v)  and 


Si  i by  performing  the  self  convolution  and  the  cross  correlation  of  each 
he;un  weighted  shape  silhouette  i t h itself  ar.d  with  all  other  i 1 lumuiatcd 
sliapes.  An  important  consideration  which  must  lie  given  ii  the  ease  of  adapt  iv 
optics  is  the  effective  differing  target  illumination  functions  which  occur  in 
adaptive  beam  forming  control  versus  the  total  illumination  which  is  sensitive 
to  the  adaptive  convergence  condition  .it  any  instant  of  time,  for  the  multi- 
dither example  considered  here  with  separate  (segmented)  adaptive  subapertures 
(a  multidither  zonal  implementation  caseB,  two  illuminating  functions  must  be 
considered.  Hiese  are  the  target  illumination  patterns  associated  with  each 
subaperture  dithered  channel,  driven  at  some  dither  frequency  <»w  for  i = 1, 

...  q dither  channels,  and  the  coherent  adaptive  control  summed  illumination 
which  results,  at  the  target  and  is  sensitive  to  the  adaptively  controlled 
instantaneous  convergence  state. 

(Considering  the  coherently  summed  illumination  pattern  first,  the 
1 suitin'.'  illumination  weighted  target  shape  silhouette  arc  checked  for  normal 
incidence  points  (target  glints  and  specular  sources).  The  self  convolution 
o!  a glint  (having  a (u,v)  plane  spatially  infinitesimal  value)  yields  i 
highly  peaked  hut  very  low  spatial  frequency  hand  contribution  to  the  PSP,  as 
hown  in  I ig.  i!  elf  convolution  of  shape  silhouettes  (being  spatially 

extended  in  the  i,  lane  is  lower  in  magnitude  but  more  extensively  spread 
in  the  frequence  In-  m ps  shown.  The  cut  off  of  this  spread  is  inversely 
related  t<  thi.  longest  dimension  of  the  (u,v)  plane  projection  of  the  target 
silhouette.  The  ci^s  correlation  of  the  dither  suhaperturc  illumination  with 
the  resultant  summed  illumination  is  in  general  spatially  broader  than  the 
sc  1 f convolved  summed  illumination  but  much  lower  in  peak  amplitude  due  to  the 
multi  dither  implementation  app roadies  which  modulates  a very'  small  fraction  o( 
the  total  beam  power  at  each  dither  frequency.  This  modulation  is  typically 
less  than  UK,  and  is  usually  in  the  order  of  2 to  3%.  With  reference  to 


the  log  scale  employed  in  Fig.  2,  the  cross-correlation  terms  will  be 
typically  15  to  17  db  below  the  glint  contributing  terns. 

As  stated  earlier,  what  is  needed  is  the  temporal  PSD.  Since  speckle 
m<  f mi  is  by  relative  target  beam  motion,  the  spatial  PSP  will  be  temporally 
modulated  by  any  target  or  beam  dynamics.  lTiis  can  he  as  complicated  as 
beam  apparent  motion  occurring  during  adaptive  control  (adaptive  search  ;ind 
lock-up  dynamics  for  a glint  point  or  a designated  target  location),  or  a 
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simply  dictated  by  target  spin  or  apparent  angular  rate  due  to  flight  dynamic- 
If  motion  effects  are  simply  associated  with  an  effective  target  rotation,  , 
then  the  temporal  PSD  can  be  simply  derived  from  the  spatial  PSP  by  multiplying 
the  self- convolution  components  by  20h.  The  cross-correlation  terms  are  modu 
lated  by  the  dither  carrier  frequencies , fj,  and  the  effective  p which 
generates  a frequency  spread  in  the  temporal  PSD  components  related  to  the 
inultidither  subaperature  contributions.  The  resulting  temporal  PSI)  is  what 
is  displayed  in  Fig.  2.  Only  a point  receiver  case  is  chosen.  For  any 
specific  system  configuration,  the  PSD  must  be  further  weighted  by  the 
spatial  transfer  function  associated  with  the  receiver  diameter.  Ibis 
will  result  in  a reduced  spectrum  content  in  the  speckle  induced  PSD  and 
modification  of  certain  frequency  component  peaks.  However,  in  order  to 
address  the  worse  case  condition  and  to  choose  a general  formulation,  a point 
receiver  is  more  appropriately  modeled  here.  In  this  modeling  case,  the 
important  terms  are  the  various  spectral  peaks  associated  with  the  glint, 
diffuse,  and  cross  correlated  dither  signal  like  noise  contributing  terms, 
fhese  are  shown  in  Fig.  2. 

Thus,  Fig.  2 illustrates  the  temporal  power  spectrum  found.  Ihc  temporal 
spectrum  may  now  be  used  to  find  the  necessary  state  models  as  in  !.q.  (13). 

The  procedure  is  to  find  ;ui  analytical  and  rational  expression  lor  the  temporal 
PSD  and  then  to  use  the  well  known  expression 

SyyM  * iHfjmlf.S^Oo)  (2d| 


relating  the  required  PSD  (Syy)  to  a linear  shaping  filter  (H(.io)l  driven  by 
unity  variance  white  noise,  ihe  realisation  oi  H(jui)  will  give  us  the  required 
di fferent ial  equat ions. 

ihe  PSD  for  each  of  the  speckle  processes  are  now  presented.  For  the 
glint,  it  is  assumed  that  the  PSD  is 
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where  2o?/6 , is  the  corn  frequency  amplutide  of  the  glint  and  6 is  the 

• • S ft 

correlation  constant  (equal  to  the  inverse  of  the  correlation  time).  The 
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correlation  time  is  assumed  large.  This  gives  a low  frequency  cut  off  for  the 
model  as  the  physics  depicts.  A similar  model  may  he  deduced  for  the  diffuse 
induced  component.  'Ill is  yields 
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The  dither  correlated  components  require  a different  model.  Iliese  models  are 
of  the  form  for  each  dither  frequency 


o . 


Ssisi(h>)  = 
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where  B j , i = 1 , ...,  q are  of  the  order  of  ft  , uv  , i = 1,  ...,  q arc  the 
dither  frequencies,  and  n - , i = I,  ...,  q arc  amplitude  scaling  parameters. 
These  may  he  obtained  as  an  approximation  to  the  actual  pST'  curves  found  from 
figure  2.  They  arc  adequate  rational  nnnroximants  to  the  actual  I's*' 

hq . (20)  and  spectral  factorization  yields  the  following  differential  equation 
for  the  models,  i.e.. 
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S.  = x,  , 
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whoK  .i  , ii  . mid  u-,  i = 1,  ...,  q are  zero  mean  white  ih  ise  with  un  ’ 
variance  and  intcorrel ated  with  each  other. 

The  coefficients  mav  he  determined  through  physical  knowledge,  exp  rime:.'  .1 
data,  or  by  on-line  use  at  adaptive  estimation  theory,  lho  use  ol  this  technique 
will  yield  the  optimal  estimator.  The  coefficients  may  he  learned  as  the  adap 
tive  optics  system  is  utilized. 

i ■> 


Thus,  equations  (21-24)  yield  the  temporally  varying  statistical  informa- 
tion about  the  speckle  process  in  a form  suitable  for  inclusion  into  the 
optimal  estimator.  The  next  section  discusses  the  remaining  models  necessary 
for  the  estimator  and  develops  the  estimator  equations  for  implementation. 


II. 


F.st imat ion  Hquat ions 


'Hie  basic  problem  that  is  to  be  addressed  is  that  of  estimating  the 
gradient  sector  9 I/&P  as  this  is  the  signal  necessary  to  adjust  the  parameters. 
This  section  develops  the  remaining  models  necessary  to  accomplish  tins  estima- 
tion as  well  as  the  estimation  of*  the  speckle  components,  and  develops  both 
the  continuous  estimator  and  the  discrete  estimator  for  this  class  of  problems. 

The  remaining  models  necessary  for  the  estimator  structure  are  temporal 
models  for  I (Pi  and  91/91'.  From  equation  (lOj  it  may  be  noted  that 


dj 

dt 


1L 

?r  ? 


dP 

dt 


:25] 


which  gives  the  necessary  dynamic  equation  for  I.  If  in  Hq.  (9)  the  parameter 
mav  be  written  as 


ek  = ckAt 


then  a limiting  process  yields  that 


dP 

dt 


- 31 


‘'k  3P, 


(2b) 


and  Eq . (32)  may  be  written  as 


dl  31 1 31 
dt  f k 3P~  3P 


(27) 


There  are  s<  veral  methods  that  may  be  used  to  model  the  dynamic  change  in  the 
gradient.  Tt  may  be  noted  that  since  the  second  partials  are  modulated  into 
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the  signal  (Eq.  (5)),  appropriate  processing  will  yield  the  matrix  3 7 1 / 9 P ‘ - 
Then 


or,  using  Eq.  1 33) 


d_  ( 31  \ 32  I dP 

dt  * 3P  * IF  dt 


d i 31  | _ ~ 3~  I / 3 1^  \ 
3t  \ 3P  ' ' k 3F  ' 3P  ) 


(28) 


(29  j 


Since  the  matrix  3? I/3P7  cannot  he  extracted  noiselessly,  process  noise  must  be 
added  to  this  equation  in  order  to  assure  the  proper  uncertainty  level  in  this 
equation.  This  yields 


d / 3 1 \ _ — 32 1 / 3 1 \ + 
Ut  ' 3P  ) Lk  3F  \ 3P  ' 


( SO ) 


where  w is  assumed  to  be  white  noise,  zero  mean  with  variance 


Elw(t)w  (t)l  = QwCt) 6 (t-T ) . 


Now,  another  dynamic  model  for  short  term  modeling  for  the  dynamic  change  in 
the  gradient  may  be 


± ( ) = . 
dt  ' 3P  I 


(31  I 


where  each  element  of  v is  a process 
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with  initial  condition 


v. 

i 


> 0 


-4 


1 


14 


where  uy  is  zero  mean  unity  variance  white  noise.  The  correlation  time  is 
chosen  on  the  order  of  the  expected  convergence  time.  The  initial  condition 
is  chosen  to  represent  the  negative  change  in  the  gradient  for  a convergent 
process.  The  choice  of  correlation  time  as  suggested  will  yield  a model  that 
gives  a decay  of  the  gradient  from  some  positive  value  to  a fluctuation  about 
zero  caused  hv  the  process  noise.  This  model  for  short  time  intervals  allows 
for  a dynamically  varying  gradient  vector.  This  effect  will  he  studied  in 
more  detail  when  the  estimator  has  been  developed. 

lqs.  (28  32J  yield  the  necessary  dynamic  models  for  the  estimator.  It 
may  be:  noted  that  this  system  of  differential  equations  may  be  placed  into  the 
form  of  hq.  118)  where  the  state  vector  x is 

x1  = U(P),  VW1,  s,  , . ..s  , x2j,  . ..x?ci>  v ,T  ,g! 1 . (88 1 

The  measurement  equation  may  he  written  in  the  form  of  hq.  (12)  as 


= k(z)U  ,1 ,1  ,. . . ,1  ,0,0, 1,1 1 + r 


where  s is  the  vector  consisting  of  the  elements  s s ar.d  x,  is  the 

q 

vector  consisting  o(  the  elements  x„  , x.  , ....  x,  . 

• 1 ' 7 •(( 

The  structure  as  given  in  hqs.  24-82  may  be  applied  to  hq.  (15)  to  yield 
the  estimator 


where  the  gain  K is  computed  via  liq.  (I”)  where 
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with  the  element  P,  being  the  covariance  of  the  estimation  error  for  i 
diagonals  P,,j  being  the  covariance  of  the  estimation  error  for  71,  am 
forth.  The  cross  terms  represent  correlations  between  th  a-  .ibles. 
row  motor  h in  il|.  ! 1 " ) is 


h = kr:)!!,:!1, 1,1,. ..,1,0,0,1,11  . 


The  matrix  P mav  be  calculated  from  Eq.  (18 J where 


0,2ekV~1 2ckVqI,  0 0,  0 

o,  o o , n o,  o 

o,  o o , o o,  n 

o,  n o , o o,  o 

0,  0 0 , 0 0,  1,  0 

o,  o , n , 0 , 0,  0,  1 

0,  0 0 , 0 0,  0 

' in,  0 0 ,-C!  ■(),..  .0,-8.  ,0 


and  ()  equal  to  a 2q  + 2 identity  matrix.  It  may  be  noted  that  jf/9x  is  a 
sparse  matrix.  This  fact  may  be  i sed  to  eliminate  many  of  the  interconnections 
and,  thus,  simplify  the  equations. 

This  structure  will  yield  the  estimate  of  the  physical  variables  to  be 
estimated.  Higher  order  approximation*  to  the  conditional  mean  may  be  obtained 
as  well  as  the  partial  integral  equation  of  the  conditional  probability 
density  function.  However,  t he  complexity  of  the  structure  is  not  worth  the 
additional  effort.  Titus,  the  above  yields  the  extended  Kalman  estimator.  The 
next  section  considers  a control  philosophy. 

IV.  Control  Philosophy 

The  problem  is  to  maximize  the  intensity  at  the  far  field  by  adjustment 
of  the  pertinent  system  parameters.  These  parameters  may  include  phase  biases 
in  segmented  mirrors  or  a deformable  mirror  and  the  focal  length  of  a tele- 
scope among  other  parameters.  The  intensity  is  a stochastic  process  duo  to 
perturbations  of  this  phenomena.  The  rerformance  index  J is  given  as 

• J(P)  = max  E{I(p)  : f 40) 


[* 
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where  ^ is  a constraint  set  for  P and  !:{•}  denotes  the  expected  value  over  all 
realizations.  Ihe  condition  for  optimization  is  that  the  stochastic  gradient 
must  he  equal  to  zero,  i.e., 
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E {VI}  = 0 


(41  ) 


where  interchange  of  expectation  and  differentiation  has  occurred.  Ihe  estima- 
tion algorithm  yields  a best  estimate  of  this  stochastic  gradient.  Thus,  the 
algorithm  to  adjust  the  parameters  in  order  to  accomplish  the  hill  climbing 
to  reach  the  stationary  point  as  depicted  in  Fq.  (9)  is 
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* P 
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where  is  a given  sequence  {r,  + j «... } to  be  chosen  and  VI  is 

the  estimator  estimate  for  V I. 

I1 

ihe  gain  sequence,  must  be  chosen  to  assume  convergence  of  the  algorithm. 
Furthermore,  it  is  desirable  that  the  convergence  by  accelerated.  Hie  method- 
ology for  choosing  the  gain  sequence  comes  from  stochastic  approximation  [9]. 

(This  is  not  an  exhaustive  bibliography.  Sec  [10J  for  a more  complete  biblio- 
graphy.) This  paper  is  not  intended  to  survey  all  existing  stochastic  approxi- 
mation methods.  However,  it  is  sufficient  to  point  out  that  the  gain  sequence 
need  not  he  exceedingly  complex.  For  example,  if  the  stochastic  gradient 
hi  1 1 P) | PI,  where  it  is  conditioned  on  P,  is  such  that 


IMF  l ( P-  P° ) 1 F.{  I ( P)  | P)  ] 0,  V*  0 (43) 

t:  * | | P-P°  | | <r " 1 

and  if  the  estimated  gradient  is  such  that 

E(!|v‘l | \7}  < hi l + ( |P-P°| !?1 , h > 0 (44) 


where  P°  is  the  value  of  the  parameters  maximizing 
sequence  of  gains  that  satisfy 


the  intensity, 


then  the 
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will  yield  aii  algorithm  convergent  in  mean  square  and  with  probability  one, 


Prob  j Jira  l>  = P°|  = 1. 
I k-*®  k 1 


The  use  of  accelerated  convergence  methods  can  be  easily  formulated  trom 
existing  control  literature.  The  philosophy  of  using  stochastic  approximation 
was  given  in  this  section.  The  next  section  considers  the  use  of  the  estimator 
for  tiie  autofocus  problem. 

V . l.st  Lna  t i oi  i in  Auto  focusing 

An  active  autofocusing  scheme  has  been  described  in  reference  | 1 l ] . 
Basically  this  scheme  consists  of  adapting  the  focal  length  by  sinusoidally 
perturbing  the  distance  between  the  secondary  mirror  and  the  primary  mirror  of 
a Cassegrain  telescope.  This  adjustment  continues  until  the  intensity  is 
maximized.  The  scheme  in  [11 | can  he  classified  as  a sinusoidal  perturbation 
adaptive  controller  as  previously  derived  in  this  paper.  This  section  gives 
the  results  of  application  of  the  estimation  to  this  problem. 

The  measurement  equation  can  be  functionally  written  as 


y - ktzMhdl  + jjj  sin  + + T + g)  + n 


where  .1  is  the  distance  between  the  primary  and  secondary  mirror  (see  Fig.  3). 
The  models  are  as  given  previously  for  the  speckle  return,  intensity,  and 

gradient . 

The  simulations  were  conducted  bv  using  ""11  known  analytical  expressions 
for  the  intensit>  of  a focused,  aperatured  Gaussian  laser.  Thus,  the  intensity 
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at  d was  calculated  by  use  of  the  analytical  expression  for  1 in  reference  1 12  | 
and  the  expression  relating  the  focal  length  to  d.  An  analytical  expression 
for  3! /3d  may  he  obtained  from 


91  _ 91  3 f 
od  9?  3d 


(.48) 


where  9I/3f  may  be  explicitly  calculated  from  reference  (12|  and  3f/9d  mav  be 
easily  determined.  The  speckle  noise  was  simulated  using  the  shaping  filter 
model  as  in  equation  (24)  with  appropriate  white  driving  noise.  These  are 
used  to  form  the  measurements  in  the  simulations.  The  estimator  is  then  used 
to  determine  the  best  estimates  for  the  gradient  and  speckle.  Die  square 
roots  of  the  diagonal  elements  of  the  covariance  matrix  gives  the  measure  of 
the  standard  deviation  of  the  estimation  error.  These  are  plotted  in  Figures 
a.  In  addition,  20(1  Monte  Carlo  runs  were  conducted  to  ascertain  the  effects 
o!  the  nonlinearities  in  the  system  dynamics.  See  Figure  4.  Figure  7 yields 
the  center  value  of  the  separation  distance  between  the  primary  and  secondary 
mirrors. 

The  modeling  approach  is  thus  based  on  Fq.  (12),  where  klz)  is  an  e f feet i v\ 
range  and  target  scattering  dependent  intensity  scaling  factor,  I(P),  is  based 
on  the  autofocus  case  discussed  in  Ref.  [12],  ;ind  speckle  induced  factors,  s^, 

T and  g are  derived  for  the  target  range  and  dynamics  associated  with  that 
example.  The  inclusion  of  the  speckle  related  interactive  effects  has  not 
been  considered  earlier,  and  allows  the  autofocus  example  to  be  more  realisth 
in  terns  of  a practical  implementation  ex amp  1 e . 

Hie  physical  variables  used  in  the  simulations  were 


zv  = 1.1  KM 

iN 

fN  = son  m 

r = 0.4  V (secondare  radius  of  curvature) 

i = 1.98  (nrimarv  radius  of  curvature) 

? ) 

d^  • 0 . "’912  1'  (mirror  spacing) 
d = 0.02  M (dither  length) 
f.  - IK  Hz  (dither  focus  rate) 

u 

. , = 2rfd  = 2U00  » 

u 

a = 0.0589  M 
b - 0.299  M 
\ = in  ft  mM 


2( 


fc.j.  =■  target  apparent  brightness  width 
= 1 M 

0^  = 10  mr/s  + 100  mr/s 

apparent  angular  rotation  rate 

B = 5 Hz 
g 

ST  = 30  Hz 

ed  = 5 112 

B = 5K  Hz 
v 

a = or  - 10 
g 1 
o,  = 100 

a =0.1. 

v 

Hie  variance  due  to  speckle  (oj)  at  dither  carriers  is  large  due  to  changes  in 
apparent  rotation  and  apparent  asoect  angle  changes  with  I light.  The  large 
magnitude  for  s indicates  Little  temporal  correlation  from  update  to  update 
for  the  correction  loop.  Ihe  estimator  update  rate  was  -IK  llz  while  the  dither 
focus  was  driven  at  f j = IK  llz. 

It  is  assumed  that  the  estimator  uses  a range-doppler  tracker  for  handovei 
giving  a value  of  range  and  a value  of  the  rotation  rate  as 

0 _ s . , x/ ■? 

donpler  IN 

where 

•TN  - U-44  x/2b) 

and  I'y  may  he  calculated  from 


Hie  problem  of  an  uncertain  rotation  rate  will  not  he  treated  here.  It. 
suffices  to  -iv  that  adaptive  estimation  techniques  mn>  be  used  to  adapt  upon 
this  variable. 


Returning  noiv  to  the  estimator  equations,  it  can  he  noted  that  in  order 
to  estimate  g,  'i  and  s,  infonnation  must  be  available  oi  the  target  apparent 
angular  motion,  9,  and  brightness  length,  Jt.j,.  The  latter  can  be  derived  from 
the  previous  consideration.  Focus  estimation  can  be  determined  from  the 
dither  focus  control  servo  (values  of  f (dQ)  which  optimize  l(.z,f).  Angular 
rates  can  be  derived  from  either  a doppler  measurement,  or  speckle  induced 
frequency  components  it  the  sensor.  1 he  autofocus  simulation  assumed  that 
initial  handoff  is  from  a range-doppler  tracker,  providing  z^(oj,  0^(0), 
f^.foj  which  are  used  to  start  the  estimator.  The  complete  description  of  the 
estimator  then  becomes 
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The  performance  was  compared  with  a simple  (nearh  ideal  phase  lock 

loop. 

VI.  Conclusions 

This  paper  generalizes  the  multidither  adaptive  optics  concepts  b> 
showing  that  they  belong  to  the  class  of  sinusoidal  perturbation  adaptive 
control  problems.  Ihis  analysis  yields  insight  into  the  higher  order  fre- 
quency terms  as  these  can  be  easily  seen  to  be  related  to  the  second  deriva- 
tive of  tile  intensity,  ibis  signal  is  corrupted  by  speckle  effects  as  well 
as  detector  noise,  lhe  speckle  components  are  modeled  in  state  space  formed 
by  the  use  < ' the  temporal  power  spectral  density  derived  in  this  paper.  The 
models  arc  then  used  to  f n an  extendea  kal.nai.  estimator  to  estimate  the 


gradient  vector.  This  separates  the  gradient  from  the  speckle  components. 

l’he  estimator  may  he  used  in  a high  bandwidth  digital  signal  processor 
similar  to  digital  radar  processors  in  order  to  estimate  the  required  quant i 
ties.  It  may  he  also  used  in  an  analog  realisation.  This  may  lead  to  an 
elimination  of  conventional  electronics  to  accomplish  the  necessary  detect ic 
while  eliminating  the  speckle  noise  from  the  measurement,  iliis  yields  a 
control  variable  that  is  less  noise  and,  thus,  better  for  feedback  purposes. 
The  control  philosophy  using  stochastic  approximation  is  briefly  discussed. 
I'll  is  yields  an  excellent  convergent  algontlim.  The  estimation  ilgoritiuus 
were  applied  to  the  autofocus  system.  Stable  convergence  occurred  in  the 
presence  of  speckle  induced  competing  signals.  A similar  case  was  tested 
with  no  estimation  employed,  leading,  in  general,  to  poor  convergence  for 
Monte  Carlo  simulated  runs. 
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Power  spectral  density  taxing  into  account  target  gum 
(specular  regions),  distributed  diffuse  regions  and  cross 
correlated  adaptive  subapertui'e  factors.  Only  a point 
receiver  is  considered. 
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Fig.  4.  Target  Intensity  - Autofocus  Performance  - Fixed  .Aspect. 


diffuse  domain  spectrum,  P., ; and  glint  domain, 
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